Criticality of Electron- Nucleus Cusp Condition to Local Effective 

Potential Energy Theories 

Xiao- Yin Pan and Viraht Sahni 

Department of Physics, Brooklyn College of the City University of New York, 
2900 Bedford Avenue, Brooklyn, New York 11210, 
and The Graduate School of the City University of New York, 
360 Fifth Avenue, New York, New York 10016. 
(Dated: February 2, 2008) 

Abstract 

Local(multiplicative) effective potential energy theories of electronic structure comprise the 
transformation of the Schrodinger equation for interacting fermi systems to model noninteract- 
ing fermi or bose systems whereby the equivalent density and energy are obtained. By employing 
the integrated form of the Kato electron-nucleus cusp condition, we prove that the effective electron 
-interaction potential energy of these model fermions or bosons is finite at a nucleus. The proof 
is general and valid for arbitrary system whether it be atomic, molecular, or solid state, and for 
arbitrary state and symmetry. This then provides justification for all prior work in the literature 
based on the assumption of finiteness of this potential energy at a nucleus. We further demonstrate 
the criticality of the electron-nucleus cusp condition to such theories by example of the Hydrogen 
molecule. We show thereby that both model system effective electron-interaction potential ener- 
gies, as determined from densities derived from accurate wave functions, will be singular at the 
nucleus unless the wave function satisfies the electron-nucleus cusp condition. 



I. INTRODUCTION 



Consider the time-independent Schrodinger equation for atoms, molecules, or solids: 

H^{ri,...rN) = E^{ri,...rN), (1) 

where \l/(ri, ...rjy) is the system wave function, E the energy, and the number of electrons. 
(We suppress the spin coordinate throughout the paper.) The electronic Hamiltonian H in 
atomic units is 

^ = E V.^ + E -(r.) + I E rr^H ' 
i i ^ i,j '■jl 

where the first and last terms represent the electron kinetic and electron interaction potential 
energy operators, respectively, and v{r) the operator representing the potential energy of 
the electrons due to the external charge Z of the nuclei at Ha'. 

-W = Err^. (3) 

The electronic density p(r) of the system is defined as 

p(r) = N J V'*(r,r2...rAr)?/'(r, r2...rAr)(ir2...(irAr. (4) 

At the coalescence of an electron with a nucleus, the external potential energy operator v{r) 
is singular. For the wave function to satisfy the Schrodinger equation and remain bounded, 
it must satisfy the Kato [l| electron nucleus cusp condition at each nucleus. The cusp 
condition is usually stated in differential form iLl. In integrated form, the statement of the 
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cusp condition for arbitrary state of the system is [2, |3|| 

V^(r, r2, ...r^) = ^^(0, r2, ...r^).(l - Zr) + r ■ a(r2, ...r^) + 0{r^), (5) 

where the vectors rj now represent the positions of the electrons from a particular nucleus. 
The vector a(r2, ...ttv) is undetermined. The derivation of the integrated form of the cusp 
condition does not involve any boundary conditions far from the nucleus. Hence, it is valid 
at the nuclei of atoms, molecules, and periodic and aperiodic solids. 

In addition to solving the Schrodinger equation(l), electronic structure is also determined 
by calculations performed within the framework of local(multiplicative) effective potential 
energy theories. Representative of such theories are Kohn-Sham density functional theory 



(KS-DFT) P], the Optimized Potential Method (0PM) and quantal density functional 
theory (Q-DFT) j^. The basic idea underlying such theories is the transformation of the 
Schrodinger equation to that of the model S system of noninteracting fermions whereby the 
equivalent density p(r) and energy E are obtained. The differential equation describing the 
model fermions is 

[-^V^ + Vs{r)]Ur) = Q0.(r); 2 = 1,2, ..N, (6) 
where Vs{r) is their effective potential energy: 

Vsir) =v{r)+Vee{r), (7) 

with fee(i') their effective electron-interaction potential energy. The potential energy fee(i') 
accounts for electron correlations due to the Pauli exclusion principle and Coulomb repulsion, 
and Correlation-Kinetic effects which are a consequence of the difference in the kinetic energy 
of the interacting and noninteracting systems. The various theories differ in their definitions 
of the potential energy fee(i')- In KS-DFT, the energy ii^ is a functional of the ground state 
density p(r), and fee(r) is defined as the functional derivative of its KS electron-interaction 
energy E^^^[p] component. In the 0PM, this energy component is considered a functional 
of the orbitals 0i(r) : E^^^ = E^^^[(f)i]. The potential energy fee(r), defined as the functional 
derivative of E^^^[(f)i], is obtained by solution of an integral equation in conjunction with 
the differential equation (6). Within Q-DFT, the potential energy fee(i") is the work done 
to move the model fermion in a conservative field J-'s{r). The density p(r) equivalent to that 
of the interacting system Eq. (4) is determined from the orbitals 0i(r) as 

p{r) = j:m^)\"- (8) 

i 

The energy is determined from the energy functionals, or in terms of the components of 
the conservative field. The highest occupied eigenvalue of the differential equation is the 
negative of the ionization energy 

An understanding of the structure and general properties of the potential energy fee(r) 
that allows for the transformation from the interacting to the noninteracting model system 
is therefore of importance. The structure is also of significance for the evaluation and 
construction of approximations within the various S system theories. A key aspect of this 
structure, and one that has been controversial, is whether the potential energy Wee(i') is 



finite or singular at a nucleus. The principal manner by which the structure is determined 
is one that employs methods jl^ which assume knowledge of the exact density p(r). The 
densities in turn are obtained from correlated or configuration-interaction wave functions 



that are 
et al 



,1C J 





lighly accurate from the total energy perspective. Work on the He atom by Smith 



Davidson , and Umrigar and Gonze jl3], show fee(r) to be finite at the 
nucleus. Almbladh and Pedroza Q|, on the other hand, showed it to be singular there. 
Additional work on light atoms show it to be either finite ^ or to diverge at the 



nucleus Il6l. For the determination of its structure for few- electron molecular systems 
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(17| such as H2 and LiH, the potential energy fee(r) is assumed finite at each nucleus. 
Expressions for fee(r) at a nucleus have also been derived [l8|, but once again they are 
msed on the assumption that it is finite there. In various approximations with KS-DFT 
the potential energy fee(r) also diverges at the nucleus. 

The controversy was resolved for closed shell atoms and open shell atoms in the central 
field approximation by Qian and Sahni [2^ who proved analytically that fee(r) is in fact 
finite at the nucleus. Furthermore, they were the first to show that this finiteness was a 
direct consequence of the electron-nucleus cusp condition. In their proof they employed 
the differential form of the cusp condition. In this form of the cusp condition, the angular 
dependence is integrated out. Hence, its application constrains their proof to spherically 
symmetric systems. 

In this paper we generalize the conclusion of Qian and Sahni and prove that fee(r) is finite 
at a nucleus independent of the type of system (atomic, molecular, or solid state), and of 
the system state and symmetry. Our proof employs instead the integrated form of the cusp 
condition, and it is for this reason that the result is valid for systems of arbitrary symmetry. 
The proof too is distinctly different. Ex post facto, the proof thus provides justification for all 
the work on the determination of f ee(r) based on the assumption that it is finite at a nucleus. 

It is also possible to transform the Schrodinger equation to that of the B system of 
noninteracting bosons such that the density p(r) and energy E of the interacting system is 



once again obtained 



In this local effective potential energy theory, the density amplitude 
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p(r) is determined directly. The differential equation for the model bosons is 

1 



[-lV^ + VB{v)]^p{v)=^^^p{v), (9) 

where ^^(r) is their effective potential energy: 

VB{v)=v{v)+vi{v), (10) 

with f£(r) the corresponding effective electron-interacting potential energy. The potential 
energy f£(r) accounts for Pauli and Coulomb correlations, and Correlation-Kinetic effects 
due to the difference in kinetic energy of the interacting fermion and noninteracting 



boson systems [2l|. Once again, within density functional theory, f£(r) is defined j8| as a 
'unctional derivative of an electron-interaction energy functional E^^[p\, whereas in Q-DFT 
2ll | it is the work done in a conservative field JFg(r). The energy is determined from the 
total energy functional or in terms of the components of the conservative field. The single 
eigenvalue /i is the chemical potential or the negative of the ionization energy. 

In this paper we also prove that the potential energy w£(r) is finite at a nucleus. The 
proof is again general and valid for systems of arbitrary state and symmetry, and also 
employs the integrated form of the cusp condition. 

As a second component to the paper, we demonstrate the criticality of the electron- 
nucleus cusp condition to the finiteness of the potential energies fee(r) and f£(r) at a 
nucleus by application to the Hydrogen molecule. It becomes evident thereby that densities 
derived from wave functions that do not satisfy the cusp condition lead to potential energies 
that are singular at a nucleus, irrespective of how accurate the wave functions may be from 
an energy standpoint. 



II. CRITICALITY OF CUSP CONDITION APPLICATION TO THE HYDRO- 
GEN MOLECULE 

For two electron systems in their ground state such as the Helium atom, Hooke's atom, 
and Hydrogen molecule, the S and B systems are equivalent. This is because the S system 
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orbital is then 0j(r) = yp(r)/2, i = 1,2. Hence, the demonstration of the significance of 
the cusp condition as apphed to the Hydrogen molecule is equally valid for both systems. 
Inverting equations (9) or (6), we then have for the Hydrogen molecule 



Ke{^)=^^+ 'r— -v{r). (11) 



(In this example, the S system differential equation has only one eigenvalue.) It is evident, 
therefore, that the singularity in v{r) at each nucleus must be cancelled by the ^ v^^^ 



term in order for Vpp(r) or f£fr) to be finite there. 



In our calculations we employ the accurate gaussian geminal wave function of Komasa 



and Thakkar 



22j | for the Hydrogen molecule which is spin free and of the form 



^ 150 
^(n, ra) = t(1 + A2)(l + Pab) E Cfc0A: (12) 

^ k=l 



in which 



2 2 2 2 2 

(j)k = exp{-ak\ri -Ra\ -Pk\ri-Rb\ -Cfc|r2-R.a| -?7fc|r2-Ra| -7fc|ri-r2| ), (13) 

where for j G {1,2} are the position vector of the electrons, Hj for j G {a,b} are 
the position vectors of the nuclei, P12 and Pab are permutation operators that interchange 
the electronic and nuclear coordinates, respectively, and Ck,ak, Pk,Ck,Vk, Ik are variationally 
determined parameters subject to the square- integrability constraint 

(ftfc + Pk) iVk + Cfc) + 7fc(afc + Pk + Vk + Ck) > (14) 

for each k. The individual exponential parameters are allowed to become negative as 
long as square-integrability is satisfied. The ground state energy obtained with this wave 
'unction is E = —1.174475313 a. u. and that of the most accurate correlated wave function 
2^ is E = —1.174475668 a.u. Thus the energy obtained by the gaussian wave function is 
accurate to the sixth decimal place. It is well known that such gaussian geminal or orbital 
wave functions do not satisfy the electron-nucleus cusp condition, and it is for this reason 
we employ this wave function. Additionally, in contrast to wave functions such as the 
Kolos-Roothan type wave function 0, 13; calculations are analytical. 
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In Fig. 1 we plot the density along the nuclear bond axis with the two nuclei on the z 
axis at R — ±0.7a.u.. It is evident that the density is very accurate right up to the nucleus, 
and on the scale of the figure appears to possess a cusp at each nuclear position. However, 
in magnifying the scale as in Fig.2, We see that there is no cusp as expected, and that the 
density is smooth across the nucleus. 

In Fig. 3 we plot v{r) and ^ ^ along the z axis about one nucleus. Whereas v{r) is 

singular as expected, the — X — term is finite at the nucleus. Therefore, the singularity 
in v{r) is not cancelled. Hence, although the wave function is very accurate from the 
perspective of the ground state energy, the fact that it does not satisfy the electron-nucleus 
cusp condition leads to the potential energy Vee{r) (or i'£(r) ) being singular at each nucleus. 
In determining these potential energies from accurate densities , it is therefore imperative 
that the densities be obtained from wave functions that satisfy the electron-nucleus cusp 
condition. 

Ill . PROOF 

We first prove that the potential energy Vee(i") is finite at a nucleus. Employing the 
integrated form of the electron- nucleus cusp condition for the wave function Eq. (5), the 
electron density near a nucleus as obtained from the definition of Eq. (4) is 

3 3 

p{r)^p{0)[{l-Zr)^ + 2{l-Zr)J2Bkrk+ E rirmAi^ + ..], (15) 

fe=l l=l,m=l 

where Bk = J akijj{0, r2, ...rjv)~^(ir2...drAr and Aim = I aiamtpiO, r2, ...rjv)~^dr2...(irAr are con- 
stants, and Tfc, ttk etc., components of the vector r and a. It follows then that 

3 3 

fe=l l=l,m=l 

= ^/^)[l-Zr + B-r + 0(r% (16) 
where in the second step we have retained only terms of 0{r). 

Inverting the S system differential equation we obtain for any occupied orbital 0i(r) the 
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FIG. 1: The electron density p{r) of the Hydrogen molecule along the nuclear bond axis in atomic 
units(a.u.). The nuclei are on this axis at ±0.7a.u.. The density is determined by the wave function 
of Eq.(12). 

expression for the electron-interaction potential energy Vee(i") as 

^ee(r) = Q + ^|^-^(r). (17) 
Next, we rewrite the orbitals 4>i{r) as 

Ur) = ^Mr), t = l,2,...N, (18) 
where the coefficients Ci(r) satisfy 

N 

Y^c,{vr = i. (19) 

i=l 

This definition of the 0i(r) is consistent with Eq. (8). Expanding the coefficient Q(r) about 
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FIG. 2: The electron density p(r) of the Hydrogen molecule near a nucleus as determined by the 
wave function of Eq.(12) in atomic units(a.u.). The nucleus is indicated by the large dot on the 
axis. 

the nucleus we obtain 

Q(r) = Q(0) + VQ(0)-r + 0(7-2) 

= Ci{Q)[l + B-r + 0{r% (20) 

where D = Vcj(0)/cj(0) is some constant vector. Inserting Eq. (16) and (20) into Eq. (18), 
we obtain the expression for the orbitals 0i(r) near the nucleus as 

0,(r) = ^/^)ciiO)[l - Zr + (B + D) • r + Oir% (21) 




y2 /~7yT 

FIG. 3: The external potential energy v{r) and the function _ , . , about a nucleus of the 



Hydrogen molecule in atomic units(a.u.). The density p{r) is determined by the wave function of 
Eq. (12). The singularity of the external potential energy at the nucleus is not cancelled by the 



function. 



Now the expression for the external potential energy near the nucleus is 

= (22) 

where the sum is over all the other nuclei. At the nucleus, the term — ^ is singular, whereas 
the other terms are constants. Prom Eq.(17) it is evident that this singularity must be 
cancelled by the \ f}y} term. 
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Consider the term ^^'^^^ near the nucleus with 0i(r) given by Eq. (21). We have 
V^r — 2/r and V^(B + D) • r 0. After acting by and taking the hmit as r — > 0, terms 

of O(r^) lead to constants while higher order terms vanish. Thus, near the nucleus, the 
term is —Z/r plus some constant, and therefore in this limit the singularity of the 

external potential energy is cancelled. Therefore VgelO) is finite. 

The proof of the finiteness of the B system potential energy ^^(r) at the nucleus is along 
the same hues as above. Substitution of Eq.(16) into the expression for ^^^(r) of Eq.(ll) 
leads to the result that f£(0) is finite. This result may also be arrived at as a special case 
of the S system proof for which Cj = 1/ \/N . 



IV. CONCLUSION 

We have proved that the effective electron-interaction potential energy of model nonin- 
teracting fermi and bose systems that reproduce the density and energy of an interacting 
fermi system in an external field is finite at a nucleus. The proof is valid for arbitrary 
state and symmetry of the interacting system. The finiteness of this potential energy at 
the nucleus in these local effective potential energy theories is a direct consequence of the 
electron-nucleus cusp condition. Since the cusp condition holds for both ground and excited 
states of the interacting system, the effective electron-interaction potential energy is finite 
at the nucleus whether the noninteracting fermions are in a ground or excited state. The 
proof is general, and does not distinguish between ground and excited states of the model 
system. The noninteracting bosons are, of course, always in their ground state. 

The study of the structure of these effective potential energies via densities derived from 
accurate wave functions now has an important proviso. These wave functions must satisfy 
the electron-nucleus cusp condition. Otherwise, the potential energies will be singular at a 
nucleus, thereby leading to erroneous conclusions regarding their structure. 
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